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Abstract. We report CO observations towards MS0735. 6+7421 a distant cooling core galaxy, and towards M87, the nearest 
cooling core in the center of the Virgo cluster. Both galaxies contain radio cavities that are thought to be responsible for the 
heating that can regulate or stop the cooling of the surrounding gas. In this feedback process, there could still be some gas 
cooling along filaments, along the borders of the radio cavities. Molecular gas is known to exist in clusters with cooling cores, 
in long and thin filaments that can be formed behind the rising bubbles inflated by the central AGN. CO emission was searched 
for at several locations along the radio lobes of those two galaxies, but only upper limits were found. These correspond to cold 
gas mass limits of a few 10 9 M o for each pointing in MS0735. 6+7421, and a few 10 6 M o in M87. This non detection means 
that either the cooling is strongly reduced by the AGN feedback or that the gas is cooling in very localized places like thin 
filaments, possibly diluted in the large beam for MS0735. 6+7421. For M87, the AGN heating appears to have stopped the 
cooling completely. 
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1. Introduction 

Giant ellipticals in galaxy cluster centers live in a complex en- 
vironment. In some cluster cores the intra-cluster medium is 
known to cool down, but how much of this gas is able to fuel the 
central galaxy is still not clear. The hot gas is no longer thought 
to be flowing directly towards the central galaxy. The regulation 
of the intra-cluster medium (ICM) cooling involves some re- 
heating mechanisms like jet feedback from a central AGN and 
it becomes observable by the formation of X-ray emitting gas 
cavities. This feedback mechanism produces mass outflows, 
weak shocks and heating ( McNamara & Nulsenl2 007). but also 
cooling, in long and thin filaments. Indeed, some very cold gas 
(10-100 K) has be en observed in CO emission in many cool- 
ing core systems fedgdl2001t ISalome & Combesll2003l) . High 
resolution imaging of CO showed the existence of molecular 
filaments in Abell 1795 dSalome & Combesl20"04l) that are con- 
fined to the borders of the radio cavitie s. These filaments e xtend 
even farther outward in NGC 1275 dSalome et al.ll2008l) . The 
cold gas formation may result from the interaction of the AGN 
with the surrounding hot gas. The filaments are thought to be 
cooling filaments that can fuel the giant elliptical with cooled 
gas. More examples of this phenomenon need to be observed 
to better understand the physics of the formation of the cold 
gas. In this work, we report a search towards two cooling core 
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systems: MS0735.6+7421 and M87 in the Virgo cluster. These 
two objects are known to be the sites of very strong feedback. 
The cooling in M0735. 6+7421 very likely occurs only in the 
very central region and M87 is a poor cooling core with almost 
no cooling, so we do not expect to find much cooled gas. 

The giant galaxy MS0735.6+7421 is at a distance of 716 
Mpc (1" correponds to 3.47 kpc). In spite of its large distance, 
its giant X-ray cavities are so large, inflated by radio lobes of 
size 550 kpc, that they are resolved by the beam of the IRAM- 
30m telescope. It is then possible to probe the cold gas dis- 
tribution and compare it with the X-ray, Ha and radio maps. 
The galaxy M S0735 .6+7421 hosts a radio source 4C+74.13. 
According to dMcNamara et al. 2005 ) the AGN activity could 
have re-heated sufficiently the ICM to stop the cooling around 
the central galaxy. Figure Q] shows that the radio lobe emis- 
sion of 4C+74.13 fills the X-ray cavities of the intra-cluster 
medium. These cavities may have displaced and compressed 
the hot gas surrounding the central cluster galaxy, as predicted 
by feedback model s of cooling flows. N egative feedback mod- 
els as described by iRevaz et al.l (12008b predict that the hot X- 
ray gas when compressed may eventually cool very efficiently 
down to very low temperatures, forming cold gas filaments vis- 
ible in CO(1-0) emission. So in this galaxy, we searched for 
molecular gas emission not only in the central region but also 
far away from the X-ray cooling core region, along the AGN 
radio lobes. 
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M87 is the most nearby X-ray bright cooling core. It lies 
at a distance of 18 Mpc. X-ray mapping around M87 reveals 
a complex structure for the cooling gas, perturbed by several 
bubbles (Young et al 2002). The X-ray-cool component de- 
tected in the radio lo bes is in blobs that are connected to the ra- 
dio buoyant bubbles dGhizzardi et alj2004l). In a deepe r XMM- 
Newton observation of M87. ISimionescu et al.l (120071) find ev- 
idence that cold, metal-rich gas is being transported outwards 
also outside the X-ray arms, possibly through bubble-induced 
mixing. Many X-ray features are better be understood if there 
are repetitive AGN outbursts. Depressions appear as the rem- 
nants of earlier outbursts, while shocks are seen as brighten- 
ings within the prominent X-ray arms. As in Perseus, shocks 
may be the most significant channel for AGN energy input into 
the cooling-flow atmospheres, a nd outbursts every 3x 10 7 yr are 
sufficient to quench the cooling dForman et alj|2005l) . 




Fig.l. A 250 X 250" image of the MS0735.6+7421 cluster 



center (iMcNamara et al.ll2005l) . one arcsec is 3.5 kpc. The 1.4 
GHz radio contours of the jets emitted by 4C+74.13 are su- 
perposed. The cavities in the X-ray emission coincide with the 
radio lobe positions (200 kpc in size at the cluster redshift). The 
circles represent the 30m telescope beam at 3mm, showing the 
7 regions we observed in CO(1-0). 



2. Observations 

The observations were made with the IRAM 30m telescope 
on Pico Veleta, Spain during two different runs in summer 
2005 and summer 2006. Weather conditions were not op- 
timal. The first run (august 2005) suffered from wind and 
significant pointing corrections. The second run (may 2006) 
was done with clouds and 10- 15 mm of precipitable water 
vapor, that means higher system temperatures. We used the 
wobbler switching mode with two 3mm and two 1.3mm re- 
ceivers operating simultaneously. The beam throw was 120". 
For MS0735.6+7421, receivers were tuned to the CO(1-0) 
and HCN(2-1) emission lines, redshifted to the velocity of the 



source. The CO(2-l) emission line is not visible from the 30m 
telescope at this redshift (z = 0.216). For M87, we observed 
simultaneously the CO(1-0) and CO(2-l) ines, at 3mm and 
1.3mm. Frequent pointings and focus checks were carried out 
using bright quasars. At 3mm, we used two 512 x 1 MHz filter- 
banks. This gives a total band of ~ 1300km/s for the CO(1-0) 
line setup. In addition, we used the two 250 x 4 MHz reso- 
lution filterbanks for the 1.3mm receivers. The beams of the 
30m telescope, at 3mm and 1.3mm, are 28" (94.795 GHz) and 
13" (218.656 GHz) respectively for MS0735.6+7421 and 23" 
(114.771 GHz) and 12" (229.537 GHz) forM87. 

The signals are expressed here in main beam brightness 
temperatures. The ratio of the antenna temperature to the main 
beam temperature is T^/T m b = B e fr/F e fr with the ratio of the 
main-beam efficiency to the forward efficiency being B e ff/F e ff= 
0.75/0.95 at 3mm and 0.52/0.91 at 1.3mm (cf IRAM-30m site 
http://www.iram.es/). 

The data were calibrated with the MIRA software and re- 
duced with the CLASS90 package. Spiky channels and bad 
scans were dropped and linear baselines were subtracted for 
each spectrum. After averaging all the spectra for each line at 
each position, the data were Hanning smoothed The results are 
summarized in Tables Q] and [2] 



3. Results 

We expected the regions where the gas can cool efficiently to 
be along the radio lobes a nd X-ray cavities, as in Abell 1795 
dSalome & Combesll2004l) . This is where we looked for a past 
or still present cooling activity with CO line observations. In 
MS0735.6+7421, we observed 7 fields to search for CO(l- 
0) emission in the brightest X-ray regions in the radio lobes. 
In M87, we observed also towards 4 offsets on the borders of 
the radio lobes. Only upper limits have been obtained in both 
galaxies. 

For MS0735.6+7421 at a redshift of z = 0.216, assuming 
a Hubble constant of 71 km/s/Mpc, and a flat universe with 
cosmological parameters Q,„ = 0.27 and Qa = 0.73, the lu- 
minosity distance of MS0735.6+7421 is D L = 1060 Mpc. The 

pc 

(119971) : 



CO lu minosity L' co [Kkm s 1 pc 2 ] is given bv lSolomon et al 



L' co = 3.25 10 7 
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We computed H2 masses by using the standard conversion 
factor Mh, = 4.6 L' co M that assumes self-gravitating 
clouds. Note that this mass could be overestimated by a fac- 
tor 5 if the gas is optica lly thin and subthermally excited 
dPownes & Solomon 1998b . We applied the same formula to 
determine the CO luminosity of M87 and assumed a distance 
of 18 Mpc (z=0.00436) to this source. The integrated CO in- 
tensity Sco [Jykm/s] is: 



I 



sco - 4.95 r mb (co)^y, 



(i) 



where r m b(CO) is the main beam brightness temperature, ob- 
tained with the CO emission line, dV is the line width and 4.95 
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the telescope S/T m b factor in Jy per K. To express upper limits, 
we smoothed the spectra to 50 km/s resolution, and assumed a 
standard linewidth of dV = 300 km/s. 

3.1. MS0735.6+7421 

The analysis of the data does not show any evidence of CO(l- 
0) in any of the 7 positions investigated. 



nosity in cooling core clusters. If we refer to this relation, we 
can evaluate the expected molecular gas mass from t he known 
Hq, lu minosity coming from the central 25 kpc ( Donahue et al. 
Il992h . This gives- 6xl0 10 M G which is slightly larger than, but 
still consistent with, the lcr upper limits we found here (offset 
[0,0]). We do not find CO line emission either far away from the 
cooling cores (other offsets), where there is no optical emission 
line. 




Offset (a,8) 


Sigma 


Ico 


Time 


T sys 


M gas 


[arcsec] 


[mK] 


[Kkm/s] 


[min] 


[K] 


[10 9 M Q ] 


(0x0) 


0.8 


0.24 


90 


137 


< 12.4 


(-22 x 0) 


0.9 


0.27 


90 


135 


< 13.9 


(-22 x 22) 


0.4 


0.12 


210 


146 


< 6.2 


(-11 x -11) 


0.8 


0.24 


66 


131 


< 12.4 


(0 x -22) 


0.4 


0.12 


174 


135 


< 6.2 


(22 x -44) 


0.3 


0.09 


573 


189 


<4.6 


(-22 x 44) 


1.9 


0.57 


31 


209 


< 29.3 1 



Table 1. Summary of CO(1-0) observations, sensitivity and 
emission upper limits for MS0735. 6+7421 . All linewidths have 
been assumed to be 300 km/s. The mass upper limits are lcr 
limits. Offsets in arcsec are from the central position J2000 RA 
07:41:44.5, DEC +74:14:40. 'Note that the last offset was ob- 
served during only 3 1 min. This explains the larger upper limit 
we obtain for that region. 



At the distance of MS0735.6+7421, the CO(1-0) spectra 
could be very broad (due to the mixing of many molecular 
clouds at di fferent velocities) an d even broader in the very cen- 
tral region (ISalome et al.ll2008l) making it very difficult to de- 
tect. Stacking the emission of all the offsets together did not 
improve much the signal to noise. This means that either the 
amount of cold gas is small, or a large velocity dispersion ex- 
ists. 

For MS0735.6+7421, the 28" field of view of the 30m tele- 
scope corresp onds to a diameter of 97 kpc. Recent observations 
of NGC 1275 ( ISalome et alj|2008l) have shown that the molec- 
ular gas, in cooling core galaxies, may lie in very thin and elon- 
gated filaments of at most ~(2 kpc x 10 kpc). If this is the case 
for MS0735. 6+7421, it means a beam filling factor of 0.4% for 
such a filament. So even with a large collection of filaments, 
we expect a significant dilution effect in the antenna beam. 

The upper limits on the amount of cold molecular gas are 
summarized in Table [T] We find on average M gas < a few 
1O 9 M inside a beam of diameter 97 kpc. In a multiphase cool- 
ing flow scenario, the mass de position rate eval uated from X- 
rays is 260 M B /yr inside 30" dGitti et alJkoOTt) . For a typical 
cooling time of 10 8 yr, this means that a steady reservoir of at 
least 10 10 M Q would exist (and even more for a star formation 
rate lower than the mass deposition rate). Our total upper limit, 
obtained by summing over all the central regions observed, is 
close to a few 1O 1O M0. Considering the uncertainties on the 
mass evaluation, this means that we cannot rule out that some 
gas is indeed co olin g but not yet visible in CO . 



Fig. 2. M87 - X-ray image of M87 by lYoung et al.1 (120021) . The 
white contours are the 6cm radio emission and the dark con- 
tours show the H a nebulosities along the edges of the radio 
lobes. The circles represent the 30m telescope beam size at 
3mm. 
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(5.5 x 11) 


1.6 


258 
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< 8.6 


CO(2-l) 


(5.5 x 11) 


2.2 


258 


1020 




CO(1-0) 


(11 xO) 


1.7 


361.5 


264 


< 9.2 


CO(2-l) 


(11 xO) 


2.5 


349.5 


465 




CO(1-0) 


(33 x -22) 


1.1 


260 


330 


<5.9 


CO(2-l) 


(33 x -22) 


3.2 


260 


821 





Edge ( 2001 ) and lSalome & Combesl ( 2003 ) have shown the 



existence of a correlation between CO emission and Ha lumi- 



Table 2. M87 - Summary of CO(1-0) and CO(2-l) observa- 
tions : sensitivity and emission upper limits. All linewidths 
have been assumed at 300km/s. The mass upper limits are given 
at lcr. 



M87 has been i dentified as a very moderate cooling core 
(IForman et alJl20071) . It is possible that the presen t AGN heat- 
ing has completely quenched the cooling flow. lYoung et al 
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(2002) compared the optical emission with Chandra X-ray im- 
ages (see Fig|2]i. Filaments are observed along the edges of 
the radio lobes, which coincide with the sites where the X-ray 
gas is cooler, as in usual cooling core clusters. The filaments 
have a total Hg lum inosity of 3xl0 40 erg.s~' in 1250 arcsec 2 , 
dSparks et alJI 19931) . 

Dynamical models of hot bubbles of gas injected by the 
active galactic nuclei (AGN), rising buoyantly into the intra- 
cluster volume, indicate that the activity cycle of the AGN 
in M87 is roughly 10 8 yr dKaiserl 12003b . The largest ra- 
dio structure is likely to be the remnant of an old bubble 
(Mat hews & Brighentil 120081) . X-ray observations show a 30 
kpc-long nearly radial filament of relatively cooler gas, aligned 
in projection with the radio lobe. 

We have observed 4 points aligned along the radio lobes 
and Ha nebulosities (see Figure|2]). Only upper limits are found 
(see Table |2j. The ce nter of the ga l axy is also not detected, 
see also NGC 4486 in lCombes eTa l. (2007). This suggests that 
the feedback heating is very likely to be efficient enough to 
suppress the cooling, otherwise we should be able to observe 
cold gas falling back towards the center. 

4. Discussion and conclusions 

We observed 7 different positions along the radio lobes of 
MS0735.6+7421 and 4 positions for M87, in CO(1-0) with 
the 30m telescope. We do not detect cold molecular gas in any 
of the observed regions. Our upper limits put significant con- 
straints on the cold gas masses. We find M gas < a few 1O 9 M 
in MS0735.6+7421 and M gas < a few 1O 5 M in M87. The ab- 
sence of cold gas could be explained if strong feedback heats 
sufficiently the intra-cluster medium and prevents any cooling. 
Nevertheless, as suggested by the recent results of hig h reso- 
lution CO mapping of NGC 1275 dSalome et al.ll2008l) . if the 
cooling happens only in very thin and elongated filaments, we 
may have missed such emission in MS0735. 6+7421. The spa- 
tial and spectral dilution of a collection of CO filaments would 
make the detection difficult in a single beam of 28" (97 kpc) at 
this distance. More sensitive observations, with higher spatial 
resolution are needed to determine whether this is the case. 
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